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1.  ELECTRON  BEAM  EXPERIMENTS 


The  procedure  of  injecting  artificial  electron  beams  in  space  in  order  to 
probe  the  ionosphere  and  magnetosphere  and  to  conduct  studies  in  basic 
plasma  physics  wa3  initiated  more  than  two  decades  ago.  In  January  1969 
Hess  et  al.  [1971]  launched  a  rocket-borne  electron  accelerator  from  Wallops 
Island.  At  an  altitude  of  a  few  hundred  kilometers  the  accelerator  ejected  a 
series  of  electron-beam  pulses  up  to  1  s  in  duration  downward  along  the 
magnetic  field  line.  The  electron  energy  ranged  up  to  9.5  keV,  and  the 
current  was  varied  up  to  a  maximum  of  490  mA.  The  highest  power  pulses 
produced  auroral  rays  in  the  upper  atmosphere  which  were  detected  on  the 
ground  by  sensitive  optical  systems.  This  experiment  demonstrated  the 
feasibility  of  propagating  electron  beams  long  distances  in  space  with 
apparently  relatively  small  alterations  resulting  from  plasma  instabilities 
or  beam  propagation  problems.  In  the  intervening  years  this  technique  of 
"active"  experiments  has  been  developed  considerably  and  has  come  to  play 
an  increasingly  important  role  in  space  physics  research.  For  a 
comprehensive  review  of  this  progress  see  the  recent  article  by  Neubert  and 
Banks  [1992]. 

There  have  been  two  main  themes  to  the  active  beam  experiments.  The 
first  is  to  serve  as  a  probe  of  naturally  occurring  fields  on  either  a  local 
{e.g. ,  the  GEOS  1  and  2  experiments  [Melzner  et  al.,  1978;  Junginger  et  al., 
1984])  or  global  (e.g.,  the  ECHO  series  of  experiments  [Winckler,  1980,  1990]) 
scale.  In  recent  years,  however,  a  second  theme  has  become  increasingly 
dominant:  determining  the  detailed  physics  of  beam-plasma  interactions. 
This  involves  a  number  of  fundamental  plasma  processes,  many  of  them  of 
a  highly  nonlinear  nature:  spacecraft  charging,  beam  dynamics,  beam 
interactions  with  neutral  gas  and  ambient  plasma,  and  the  generation  of 
electromagnetic  waves  from  continuous  and  pulsed  beams. 

Active  experiments  in  space  offer  both  advantages  and  disadvantages 
compared  to  laboratory  experiments.  Among  the  advantages  are  that  the 
spatial  scales  are  much  larger  and  that  one  can  conduct  studies  free  from 
the  contaminating  effects  of  nearby  chamber  walls.  The  major 
disadvantage  is  that  diagnostic  capabilities  are  much  more  limited.  This 
latter  drawback  has  been  partially  overcome  in  recent  years  by  the 
development  of  increasingly  realistic  simulation  models  which  can  provide 
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a  global  perspective  on  the  various  processes  involved  in  the  beam-injection 
experiments.  These  models  can  thus  serve  as  an  important  complement  to 
the  actual  observations,  which  are  obtained  on  one,  or  at  most  a  few, 
spacecraft. 

In  this  report  we  review  our  progress  in  employing  three-dimensional 
electrostatic  and  two-dimensional  electromagnetic  particle  simulations  to 
investigate  the  basic  physical  processes  that  occur  during  active  electron 
beam  experiments  in  the  ionosphere. 

2.  ISOLATED-SYSTEM  SIMULATION  MODEL  IN  2D  AND  3D 

In  early  simulation  models  of  beam-plasma  interactions  [see,  e.g. , 
Dawson  and  Shanny,  1968]  it  was  customary  to  assume  periodic  boundary 
conditions  along  the  direction  of  beam  propagation.  This  meant  that  one 
was  considering  only  a  small  portion  of  a  much  longer  (strictly,  infinite 
length)  beam.  Such  a  model  is  clearly  inappropriate  for  space  experiments 
where  a  paramount  question  concerns  the  ability  of  the  beam  to  propagate 
away  from  the  near  environment  of  the  electron  gun  source.  In  simulations 
with  one  spatial  dimension,  it  is  straightforward  to  implement  a 
nonperiodic  geometry  [see,  e.g. ,  Katz  et  al.,  1986;  Okuda  et  al.,  1987; 
W inglee  and  Pritchett,  1987].  Such  models,  however,  are  restricted  to 
studies  of  field-aligned  injection.  In  this  section  we  describe  the 
development  of  an  isolated-system  electrostatic  simulation  model  in  two 
(BIG)  and  three  (BIG3D)  spatial  dimensions  [Pritchett  and  Winglee,  1987; 
Winglee  and  Pritchett,  1988;  Pritchett ,  19911. 

The  basic  configuration  of  the  model  is  illustrated  in  Figure  1  for  the 
2D  case.  The  spacecraft,  beam,  and  ambient  plasma  are  treated  as  an 
isolated  system.  The  ambient  magnetic  field  is  assumed  to  be  uniform  and 
points  in  the  x  direction.  The  spacecraft  has  an  idealized  rectangular  shape 
and  is  immersed  in  a  uniform  ambient  plasma.  It  is  located  well  away 
from  any  boundaries,  so  that  it  can  draw  return  currents  from  all  sides. 
The  solution  of  Poisson's  equation  for  the  electrostatic  potential  is  obtained 
by  convolving  the  Fourier  transform  of  the  charge  density  on  the  2D  or  3D 
grid  with  that  of  the  Green's  function  for  a  single  particle.  The  surface  of 
the  spacecraft  is  treated  as  an  equipotential,  which  is  accomplished 
through  the  use  of  the  capaticance  matrix  [Hockney,  1968].  The  beam 
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Figure  1.  Configuration  used  in  the  isolated-system  simulation  model 
in  two  and  three  dimensions.  The  rectangular-shaped 
spacecraft  is  immersed  in  a  uniform  ambient  plasma,  and 
the  uniform  magnetic  field  points  in  the  x  direction.  The 
beam  electrons  are  emitted  from  one  side  of  the  spacecraft. 
For  a  typical  rocket  experiment  the  grid  spacing  A 
corresponds  to  about  20  cm. 
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electrons  are  injected  from  one  side  of  the  spacecraft  with  a  monoenergetic 
velocity  distribution  at  any  angle  relative  to  the  magnetic  field.  For  each 
electron  injected  into  the  plasma,  the  spacecraft  is  given  an  opposite 
charge.  If  a  plasma  particle  or  a  returning  beam  electron  strikes  the 
spacecraft,  it  is  removed  from  the  system  and  its  charge  is  added  to  that  on 
the  spacecraft. 

It  is  not  possible  to  model  precisely  the  various  length  scales  and 
density  ratios  present  in  the  injection  experiments.  Thus,  the  beam 
typically  has  an  initial  width  of  a  few  centimeters,  while  the  rocket  has  a 
length  of  several  meters.  The  beam  density  at  injection  is  of  the  order  of 
107  cm-3,  while  the  ambient  plasma  density  is  only  of  the  order  of  105  cm'3. 
In  the  3D  model  the  rocket  length  is  usually  16A,  so  that  the  grid  spacing 
A  ~  20  cm.  The  minimum  beam  area  that  can  be  resolved  is  2A  x  2A  =  4A2. 
This  is  considerably  larger  than  the  actual  area.  We  shall  thus  require  only 
that  the  total  injected  current  in  the  simulation  be  the  same  as  the 
experimental  value  Io-  As  a  consequence,  the  injection  density  in  the 
simulation  will  be  smaller  than  the  actual  value.  By  matching  to  the  beam 
energy  Eb  and  Io,  we  find  that  the  dimensionless  beam  velocity  v?  s  vb/AtOpb 
in  the  simulations  is  given  by 

vb  =  (4  me  eovbfyelo)172 

=  1.15  (Eb/keV)3/4  /  do/Arnp)^2.  (1) 

Thus,  for  a  1-keV,  50-mA  beam,  which  was  typical  for  the  CHARGE  2 
[Myers  et  al.,  1989]  and  Spacelab  2  [Gurnett  et  al.,  1986]  experiments, 
vb  =  5.15.  In  contrast,  for  the  ECHO  7  experiments  [Winckler  et  al.,  1989], 
where  a  much  higher  beam  energy  was  employed  (40  keV,  250  mA),  vb  = 
36.7.  For  the  CHARGE  2B  mission  involving  a  3-keV,  2-A  beam 
[J.  Ernstmeyer,  private  communication,  1990],  vb  =1.86. 

The  typical  simulation  grid  dimensions  in  the  3D  model  are  Nx  x  Ny  x 
Nz  =  512  x  32  x  32.  The  grid  spacing  in  the  x  direction  can  be  stretched  to 
Ax  =  2A  so  that  the  simulation  can  be  run  for  longer  times  before  the  beam 
reaches  the  boundary.  This  feature  is  employed  only  when  it  is  not 
necessary  to  resolve  spatial  dimensions  on  the  order  of  the  beam  Debye 
length  X,Db  =  VbA  near  the  spacecraft. 
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The  electrons  and  ions  comprising  the  ambient  plasma  are 
represented  by  superparticles.  For  the  plasma  electrons,  the  charge-to- 
mass  ratio  is  the  same  as  for  the  beam  electrons,  and  the  charge  is  adjusted 
so  that  the  initial  density  of  the  plasma  relative  to  the  injected  beam  density 
has  the  desired  value  neo/nbo-  The  mass  of  the  plasma  ions  is  taken  to  be  100 
times  the  mass  of  the  plasma  electrons.  The  initial  velocity  distributions  of 
the  plasma  electrons  and  ions  are  Maxwellians  with  Te  =  Ti  and  a  typical 
thermal  velocity  vre  =  0.1  vt,.  With  a  superparticle  density  of  one  per  A3  cell, 
the  ambient  plasma  is  represented  by  just  over  1  million  electrons  and  ions 
each. 

3.  CHARACTERISTIC  INJECTION  TIME  SCALES 

There  are  two  characteristic  time  scales  which  determine  the  ability  of  the 
beam  to  escape  from  the  near  environment  of  its  source  [Winglee  and 
Pritchett,  1987;  Pritchett  and  Winglee,  1987J.  The  first  is  the  beam 
stagnation  time  ts.  This  is  defined  to  be  the  time  required  for  the  first  beam 
electron  to  be  brought  to  rest  under  the  influence  of  the  positively-charged 
spacecraft  for  injection  into  vacuum.  The  corresponding  distance  away 
from  the  injection  point  where  the  first  electron  comes  to  rest  is  called  the 
stagnation  distance  xs.  In  one  dimension  (the  case  of  an  infinitely  wide 
beam),  ts  =  2  ©pb'1  independent  of  the  beam  velocity,  and  xs  =  0.5  X.Db-  In  two 
dimensions,  the  stagnation  time  depends  on  vb  approximately  as  vb  1/2  and 
also  varies  somewhat  with  beam  width  and  spacecraft  size  [Winglee  and 
Pritchett,  1988].  A  typical  value  was  ts  =  9  ©pb"1,  with  xs  =  1.3  Xj)b- 

The  variation  of  ts  as  a  function  of  the  dimensionless  velocity  vb  in 
three  dimensions  is  shown  in  Table  1.  In  these  simulations  the  spacecraft 
dimensions  are  L*  x  Ly  x  Lz  =  4A  x  16A  x  8A,  and  the  ambient  magnetic  field 
strength  is  assumed  to  be  such  that  Qe/®pb  =  1.  The  most  dramatic  feature 
of  the  results  in  Table  1  is  that  the  values  for  t8  are  much  larger  than  the 
corresponding  results  in  two  dimensions.  Thus,  for  vb  =  5,  ©pbts  =  64, 
compared  to  8.3  for  a  4A  x  16A  spacecraft  in  two  dimensions.  The  reason  for 
this  substantial  change  is  that  in  three  dimensions  the  electric  field  for  a 
point  charge  falls  off  as  1/r2,  whereas  in  two  dimensions  the  decrease  is 
only  as  1/r.  It  is  apparent  that  the  stagnation  time  increases  linearly  with 
vb  .  If  we  extrapolate  to  the  value  Vb  =  36.7  corresponding  to  the  peak  ECHO 
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7  beam  energy,  we  find  (Opbtg  =  460.  The  stagnation  distance  x8  relative  to  the 
beam  Debye  length  increases  only  slightly  for  the  range  of  Vb  values  shown 
in  Table  1,  from  Xs/^Db  =  3.0  for  vb  =  2  to  3.9  for  Vb  =  8. 


TABLE  1.  Stagnation  Time  ts  as  a  Function  of  the  Dimensionless 

Beam  Velocity  vb 


Vb 

©pbts 

2 

22 

3 

37 

4 

51 

5 

64 

8 

100 

The  values  of  the  stagnation  time  depend  sensitively  on  the  true  three- 
dimensional  variation  of  the  electrostatic  field.  In  a  pseudo-three¬ 
dimensional  model  in  which  periodic  boundary  conditions  are  assumed  in 
the  y  and  z  directions,  there  will  be  repeated  images  of  the  positively 
charged  spacecraft.  These  image  charge  distributions  will  modify  the 
spatial  dependence  of  the  electric  field.  Thus,  with  a  spacecraft  length 
Ly  =  16A  and  a  periodicity  length  of  32A  in  y,  the  value  of  oopbts  for  vb  =  4  is 
reduced  from  51  to  25. 

The  second  characteristic  time  scale  associated  with  the  beam 
injection  is  the  plasma  response  time  trp.  This  is  the  time  required  for  the 
ambient  plasma  to  respond  to  the  fields  produced  by  the  injection  of  the 
electron  beam,  trp  is  determined  by  the  ambient  density,  trp  =  V3/cope,  where 
©pe  is  the  ambient  plasma  frequency.  The  relative  value  of  t8  and  t|-p 
determines  two  characteristic  injection  regimes.  When  ts  is  the  shorter 
time,  there  is  not  sufficient  time  for  return  currents  to  flow  and  neutralize 
the  spacecraft  before  the  stagnation  region  forms.  In  this  case  there  will  be 
strong  vehicle  charging  and  low  injection  efficiency.  If,  on  the  other  hand, 
trp  is  the  shorter  time,  then  the  spacecraft  will  be  neutralized  before  the 
stagnation  region  has  a  chance  to  form.  In  this  case  the  injection  efficiency 
will  be  high. 
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4.  SPACECRAFT  CHARGING 


Figure  2  shows  the  spacecraft  potential  Oc  as  a  function  of  time  for  the 
case  of  vb  =  4,  fV^pb  =  0.5,  and  for  4A  x  16A  x  4A  and  2A  x  8A  x  2A  spacecraft 
situated  in  vacuum.  With  no  ambient  return  current,  ed>c  increases 
basically  linearly  in  time  up  past  the  beam  energy.  For  the  smaller  craft  the 
stagnation  time  is  shorter  (27  CDpb'1  as  opposed  to  50  ©pb'1),  and  so  the 
potential  rises  more  rapidly.  The  charging  above  the  beam  voltage  after 
stagnation  is  greater  for  the  smaller  spacecraft,  where  the  final  spacecraft 
potential  is  some  50%  above  the  beam  value.  This  occurs  because  beam 
electrons  are  scattered  in  the  stagnation  region  onto  adjacent  field  lines. 
With  a  small  craft,  a  larger  fraction  of  these  field  lines  do  not  intersect  the 
craft,  and  so  these  electrons  will  not  be  returned  to  neutralize  the  source. 
This  results  in  a  higher  potential.  Such  charging  up  to  and  above  the  beam 
voltage  has  been  observed  during  the  SEPAC  experiments  [Obayashi  et  al., 
1985;  Sasaki  et  al.,  1986]  when  the  conductive  surfaces  of  the  orbiter  were 
out  of  contact  with  the  ambient  plasma  and  during  rocket  experiments  at 
times  of  extremely  low  ambient  plasma  density.  In  the  MAIMIK 
experiment  [Maehlum  et  al.,  19883,  the  electron  gun  was  located  on  the 
daughter  which  had  a  surface  area  of  only  1.1  m2.  This  small  area, 
together  with  the  reduced  plasma  density  of  104  cm*3,  is  evidently 
responsible  for  the  observed  charging  of  up  to  1.6  times  the  beam  voltage. 

Figure  3a  shows  the  time  history  of  the  spacecraft  potential  for 
injection  into  an  ambient  plasma  for  three  different  values  of  the  ratio 
neo/nbo  of  ambient  to  beam  density.  For  neo/nbo  =  0.25,  e<X>c  never  exceeds  5% 
of  the  beam  energy;  for  neo/nbo  =  0.10,  eOc  initially  grows  rapidly  to  10%  of 
the  beam  energy  but  then  increases  more  slowly  with  time  as  the  return 
current  flows  partially  neutralize  the  charge  accumulation  on  the 
spacecraft.  For  the  low  density  case  of  neo/nbo  =  0.025,  eOc  maintains  a 
rapid  increase  throughout  the  run  and  exceeds  the  beam  energy  for 
copbt  >  115. 

Figure  36  shows  the  time  history  of  the  return  current  to  the  spacecraft 
for  the  highest-  and  lowest-density  cases.  Since  no  beam  electrons  return  to 
the  spacecraft,  the  return  current  is  provided  entirely  by  the  ambient 
plasma.  The  total  thermal  current  to  the  spacecraft  It,  divided  by  the 
injection  current  Iq  is  found  to  be  It/I o  =  2.9  iie/nbO-  Thus  for  all  three 
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Figure  2.  Time  history  of  the  spacecraft  potential  <£c  (normalized  to 
the  beam  energy  mevb2/2)  for  injection  into  vacuum  for  two 
different  size  craft. 
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Figure  3.  Time  histories  of  (a)  the  spacecraft  potential  3>c  (normalized 
to  the  beam  energy  meVb2/2)  and  (b)  the  return  current  to  the 
spacecraft  Ir  (normalized  to  the  injection  current  Io)  for 
different  values  of  the  ambient  to  beam  density  ratio  neO  / 
nbO- 
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density  ratios  considered  here,  It/Io  <  1-  For  neo/nbo  =  0.25,  Io  is  some  40% 
larger  than  It,  and  this  small  excess  of  injection  current  can  be  balanced 
by  a  very  low-level  charging  of  the  spacecraft,  as  can  be  seen  from  Figure  3. 
For  ne()/nbO  =  0.10,  the  observed  return  current  is  about  3It,  and  this 
current  can  be  collected  as  a  result  of  a  moderate  charging  of  the 
spacecraft.  For  neo/nbo  =  0.025,  the  actual  return  current  can  only  approach 
Io  by  charging  the  spacecraft  up  to  the  beam  potential.  It  is  thus  clear  that 
the  thermal  current  It  should  not  provide  the  limit  to  the  injection  current 
that  can  propagate  away  from  the  vicinity  of  the  source.  The  degree  of 
charging  of  the  spacecraft  will  be  determined  by  the  need  to  balance  Io- 
Clearly,  the  larger  the  conducting  area  of  the  spacecraft,  the  smaller  will  be 
the  potential. 

There  is,  however,  a  limit  on  the  amount  of  current  that  can  be 
injected.  As  Io  is  increased,  the  parameter  vb  decreases,  and  hence  the 
stagnation  time  decreases.  If  t8  becomes  comparable  to  the  plasma 
response  time,  then  the  beam  will  be  brought  to  rest  before  the  return 
current  flows  can  control  the  spacecraft  charging.  The  precise  condition  is 
t8  =  2  trp  [. Pritchett ,  1991],  which  occurs  in  the  present  simulations  for  vb  «  2. 

5.  BEAM  PROPAGATION 

We  now  consider  the  propagation  of  the  beam  away  from  the 
spacecraft.  Figure  4  shows  vx-x  phase  spaces  at  late  times  for  the  case  of 
parallel  injection  for  Vb  =  4,  Qe/®pb  =  0.5,  and  a  4A  x  16A  x4A  spacecraft  for 
density  ratios  neo/nbo  =  (a)  0.25,  ( b )  0.10,  and  (c)  0.025.  The  ratio  neo/nbo  =  0.25 
corresponds  to  an  ambient  density  of  1.3  x  105  cm'3,  which  is  a  typical 
ionospheric  value.  For  the  two  higher  ambient  densities  (Figures  4a  and 
46),  the  beam  initially  propagates  away  at  the  injection  velocity.  After  a 
distance  of  ~  80A,  a  coherent  wave  structure  begins  to  appear  which  is  then 
further  amplified.  This  leads  to  the  formation  of  large  trapping  vortices  in 
which  some  of  the  beam  electrons  have  their  parallel  velocity  reduced  to 
zero.  In  the  region  of  the  first  complete  vortex  the  beam  suffers  an  abrupt 
transverse  disruption. 

Figure  5a  shows  the  potential  F  as  a  function  of  x  for  y  =  Ly/2,  z  =  Lz/2 
for  the  n«0/nb0  =0.25  case.  For  a  distance  of  -60A  in  front  of  the  spacecraft 
(the  injection  point  of  the  beam  is  x  =  256)  the  potential  is  very  small  and 
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Figure  4.  The  vx  —  x  phase  space  for  the  beam  electrons  for  injection 
into  an  ambient  plasma  with  density  given  by  neo/nbo  =  (a) 
0.25,  (b),  0.10,  (c)  0.025.  In  Figure  4a  the  time  copbt  =  130, 
while  in  Figures  4b  and  4c,  Opbt  =  140. 
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nearly  constant.  This  indicates  that  initially  the  beam  is  well  neutralized.  A 
rapid  modulation  of  4>(x)  then  grows  up  and  decays  away  over  the  range 
350  A  ^  650.  This  is  the  region  of  the  strong  trapping  vortices  (cf.  Figure 
4a).  The  peak  potentials  reach  10-15%  of  the  beam  voltage,  which  is  much 
higher  than  the  spacecraft  potential  of  3-4%.  Figure  56  shows  the 
corresponding  power  spectrum  of  the  potential  computed  at  the  point 
x/A  =  406,  y  =  Ly/2,  z  =  IV2.  The  high-frequency  oscillations  associated  with 
the  coherent  structures  peak  at  a  frequency  CD  =  0.38  CDpb  =  0.76  tope-  With  a 
wavelength  of  49A,  the  dispersion  relation  for  the  oscillations  is  found  to  be 
(D  =  0.74  kvb.  The  corresponding  relation  in  an  infinite  plasma  with  equal 
beam  and  plasma  densities  (which  is  roughly  the  situation  in  the  coherent 
region  in  the  simulations)  is  co  =  kvb/2.  The  ions  are  too  massive  to  respond 
to  these  high-frequency  plasma  oscillations.  Their  gyroradius  is  large 
compared  to  the  beam  width,  and  they  are  able  to  flow  across  the  magnetic 
field  to  neutralize  local  charge  accumulations.  This  is  the  origin  of  the  low- 
frequency  peak  near  the  ion  plasma  frequency  in  Figure  56. 

We  now  consider  the  case  where  the  beam  is  injected  at  an  angle  of 
45*  to  the  magnetic  field.  Figure  6a  shows  a  x-y  scatter  plot  of  the  beam 
electrons  at  time  ©pbt  =  100.  It  is  apparent  that  after  about  2  or  3  gyroperiods 
the  helical  structure  of  the  beam  has  been  destroyed  by  the  space-charge 
oscillations  [Winglee  and  Pritchett,  1988;  Pritchett,  1990].  This  loss  of 
parallel  coherence  is  associated  with  the  development  of  a  large  spread  in 
the  parallel  beam  velocity  caused  by  the  space-charge  fields  and  occurs  over 
a  time  scale  of  the  order  of  30  CDpb*1.  There  is  also  a  loss  of  transverse 
coherence.  Figures  66  and  6c  show  y-z  scatter  plots  for  the  beam  electrons 
in  the  range  128  <  x/A  <  188  and  248  <  x/A  <  308,  respectively.  Near  the 
source  the  trajectories  of  the  beam  electrons  form  a  hollow  cylinder.  The 
radius  of  the  cylinder  is  about  equal  to  the  beam  gyroradius  and  has  a  width 
of  about  3A,  which  is  approximately  the  beam  Debye  length.  At  greater 
distance  from  the  source  (Figure  6c)  the  cylinder  has  been  filled  in,  due  to 
the  transverse  space-charge  oscillations  excited  by  the  beam.  The  prediction 
that  the  beam  electrons  should  initially  lie  on  a  hollow  cylinder  for  the  case 
of  low  spacecraft  charging  is  consistent  with  observations  during  the  STS  3 
beam  experiments  [Banks  and  Raitt,  19881. 

Figure  7  shows  tiome  histories  of  the  elec  trie  field  Ey  and  the 
associated  power  spectra,  both  inside  the  beam  cyliner  (at  y/A  =  33,  Figures 
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Figure  5.  (a)  The  electrostatic  potential  <I>(x)  (normalized  to  the  beam 
energy  mevb2/2)  at  y  =  Ly/2,  z  =  Lz/2  at  time  copbt  =  130  for 
injection  into  an  ambient  plasma  with  density  given  by 
ne(/nbo  =  0-25.  (b)  The  corresponding  power  spectrum  of  the 
potential  computed  at  x/A  =  406,  y  =  Ly/2,  z  =  Lz/2. 
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Figure  6.  Scatter  plots  of  the  location  of  the  beam  electrons  at  time 
Cflpbt  =  100  for  injection  at  45”  to  the  magnetic  field.  (a)  x-y 
scatter  plot  for  all  electrons;  (b)  y-z  scatter  plot  for  electrons 
in  the  range  128  <  x/A  <  188;  (c)  y-z  scatter  plot  for  electrons 
in  the  range  248  <  x/A  <  308. 


7a,  lb)  and  outside  (at  y/A  =  45,  Figures  7c,  Id)  for  x/A  =  218  and  z  =  Lz/2. 
Inside  the  beam,  Ey  shows  an  increased  oscillation  after  the  beam  front  has 
passed  at  ©pbt  =  30.  The  spectrum  is  dominated  by  a  peak  slightly  below  the 
upper  hybrid  frequency  (DUH  =  (©pe2  +  f2e2)^2  =  1.4  ©pe.  It  is  presumably 
these  oscillations  that  produce  the  transverse  spreading  of  the  beam. 
Outside  the  beam,  the  intensity  of  these  high-frequency  oscillations  is 
reduced  by  about  a  factor  of  20.  Here  the  most  intense  modes  are  at  lower 
frequencies  near  the  ion  plasma  frequency. 

Thus  we  conclude  that  in  all  the  cases  considered  here  the  beam  is 
subject  to  strong  distortion  or  disruption  as  it  propagates  away  from  the 
spacecraft.  These  effects  are  expected  to  be  substantially  diminished  at  low 
altitudes  where  the  neutral  density  is  high  enough  so  that  the  effects  of 
ionization  need  to  be  included  [Winglee,  1990]. 
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Figure  7.  Time  histories  of  the  perpendicular  electric  field  Ey  (left- 
hand  side)  and  the  corresponding  power  spectra  (right 
hand  side)  for  the  same  conditions  as  in  Figure  6.  The  fields 
are  measured  at  x/A  =  218,  z  =  Lz/2  and  (a,b)  y/A  =  33  (inside 
the  beam  cylinder),  (c,d)  y/A  =  45  (outside  the  beam 
cylinder). 
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6.  WHISTLER  MODE  RADIATION  FROM  DC  ELECTRON  BEAMS 


The  generation  mechanism  of  the  whistler  waves  observed  during  the 
Spacelab  2  mission  was  not  immediately  evident.  As  noted  by  Farrell  et  al. 
[1988],  the  observed  wave  intensity,  which  was  estimated  to  be  about  3  x  10*5 
of  the  beam  power,  was  of  the  order  of  107  times  greater  than  expected  on 
the  basis  of  incoherent,  single-particle  radiation  [see,  e.g.,  Mansfield,  1967]. 
Thus  some  coherent  radiation  mechanism  must  have  been  involved.  Since 
the  beam  width  was  less  than  one  whistler  wavelength,  however,  it  is 
difficult  to  imagine  that  a  conventional  whistler  instability  could  produce 
sufficient  amplification  within  the  narrow  confines  of  the  beam.  Farrell  et 
al.  [1988]  suggested  that  the  mechanism  could  involve  coherent  Cerenkov 
radiation  from  electron  bunching  in  space  at  the  plasma  frequency 
produced  by  an  electrostatic  beam-plasma  instability. 

In  order  to  simulate  the  whistler  generation  process,  it  is  necessary  to 
use  an  electromagnetic  model  which  includes  self-consistent  magnetic  and 
inductive  electric  fields.  We  describe  in  this  section  the  results  of  two- 
dimensional  particle  simulations  which  solve  the  full  set  of  Maxwell's 
equations.  The  initial  results,  presented  in  subsection  A,  use  a  pure  initial- 
value  model  for  the  beam-plasma  interaction  [Pritchett  et  al.,  1989]. 
Subsequently,  in  subsection  B,  we  use  the  electromagnetic  version  of  the 
isolated-system  model. 

6.1.  Initial-Value  Model 

Figure  8  gives  a  schematic  diagram  of  the  initial-value  simulation 
model.  A  dilute  beam  of  field-aligned  electrons  is  assumed  to  exist  initially 
throughout  the  central  region  of  the  simulation  box.  In  the  results 
presented  here,  the  beam  density  nbo  is  taken  to  be  one  quarter  of  the 
ambient  density  neO,  the  beam  energy  is  5.7  keV,  and  the  width  of  the  beam 
is  14  m.  The  current  carried  by  the  electron  beam  must  be  balanced  by  a 
return  current  in  the  ambient  plasma.  To  represent  this  current,  we  assign 
an  initial  drift  of  -vb/3  to  the  ambient  electrons  located  in  the  beam  region 
(whose  density  is  taken  to  be  3neo/4).  Outside  the  beam  region  is  a  stationary 
ambient  plasma  of  density  neo-  The  magnetic  field  is  oriented  in  the  y 
direction  and  has  a  strength  such  that  De/ope  =  0.3.  The  simulation  system 
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Figure  8.  A  schematic  of  the  initial-value  simulation  model.  The 
initial  beam  and  return  currents  are  field  aligned  and 
confined  to  the  central  portion  (in  x)  of  the  simulation  box. 
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is  periodic  in  y  but  bonded  in  x  and  represents  a  domain  of  220  m  x  110  m. 
The  ions  are  treated  as  a  fixed  neutralizing  background. 

As  is  well  known,  the  interaction  of  a  weak  beam  with  an  ambient 
plasma  leads  to  the  excitation  of  electrostatic  Langmuir  oscillations  with 
wavenumber  k  given  by  k  «  cope  /vb-  In  the  present  simulation  this  most 
unstable  mode  corresponds  to  mode  7  in  the  y  direction.  Figure  9a  shows 
the  time  history  of  the  parallel  (y)  component  of  the  electrostatic  energy  for 
this  mode.  The  exponential  growth  is  already  apparent  by  copet  -  10,  and 
saturation  occurs  by  (0pet  -  25.  The  linear  growth  rate  is  0.23  cope'1-  The 
growth  of  the  Langmuir  oscillations  leads  to  the  well-known  formation  of 
vortices  in  the  vy-y  phase  space.  The  Langmuir  oscillations  are  confined  to 
the  finite  beam  region,  in  agreement  with  Spacelab  2  observations  [Gurnett 
et  al.,  19861.  At  saturation  the  total  electrostatic  field  energy  amounts  to  14% 
of  the  initial  beam  energy.  After  saturation,  the  Langmuir  waves  collapse 
as  the  beam-plasma  system  is  thermalized. 

Figures  9b  and  9c  show  the  time  history  of  the  z  component  of 
the  magnetic  field  energy  for  the  two  longest  wavelength  modes  in  the 
system  (modes  1  and  2  in  the  y  direction).  There  is  a  clear  exponential 
growth  phase  in  both  modes  that  becomes  apparent  above  the  background 
by  topet  -  20.  The  rapid  growth  continues  until  copet  -  30;  subsequently,  the 
field  energy  remains  at  an  enhanced  level.  In  contrast,  the  magnetic  field 
energies  for  modes  3-7  collapse  back  to  the  thermal  level  as  the  electrostatic 
waves  decay.  It  is  thus  clear  that  persistent  long-wavelength 
electromagnetic  modes  are  being  excited  on  a  time  scale  comparable  to  that 
of  the  beam-plasma  interaction. 

To  identify  the  nature  of  these  electromagnetic  modes,  we  examine  the 
power  spectra  curves  given  in  Figure  10.  These  spectra  are  measured  well 
outside  the  beam  region.  The  top  two  curves  are  for  the  y  and  z  components 
of  the  magnetic  field,  while  the  bottom  one  is  for  the  z  component  of  the 
electric  field.  In  all  three  cases  the  dashed  curve  represents  the  thermal 
level  in  the  simulation  in  the  absence  of  an  electron  beam,  while  the  solid 
curve  is  the  result  with  the  electron  beam.  The  magnetic  plots  show 
enhancements  by  factors  of  80  and  20,  respectively,  above  the  thermal  level 
in  the  frequency  range  below  Qe.  Otherwise,  there  is  no  significant  wave 
excitation  above  the  background  level.  The  long-wavelength  modes  excited 
by  the  electron  beam  are  thus  clearly  whistler  waves. 
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Figure  9.  Time  histories  of  (a)  the  y  component  of  the  electrostatic 
field  energy  for  mode  7  and  (b,c)  the  z  component  of  the 
magnetic  field  energy  for  modes  1  and  2. 
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The  wave  vectors  for  the  four  most  intense  modes  excited  in  the 
simulation  make  angles  of  27*,  45*,  27*,  and  37*  with  respect  to  the  magnetic 
field  and  have  frequencies  oi/Qe  =  0.53,  0.50,  0.72,  and  0.67,  respectively.  The 
corresponding  resonance  cone  angles  are  57*,  59*,  43*,  and  47*.  Thus  the 
most  intense  whistler  waves  are  propagating  considerably  away  from  the 
resonance  cone.  The  total  energy  radiated  as  whistler  waves  is  estimated  to 
be  0.1%  of  the  initial  beam  energy. 

The  present  simulations  show  that  the  beam-plasma  interaction 
produces  a  complex  current  structure  with  significant  Jx,  Jy,  and  Jz 
components.  While  the  largest  contribution  is  the  mode  7  component 
associated  with  the  dominant  Langmuir  wave,  there  are  also  substantial 
contributions  in  modes  1-6.  The  transverse  components  Jx  and  Jz  arise 
from  the  existence  of  the  electrostatic  field  in  the  x  direction  which 
accelerates  the  electrons  perpendicular  to  the  magnetic  field  and  produces 
E  x  B  drifts.  Figure  11  shows  time  history  plots  of  IJy(ky)l2  and  IJz(ky)l2  for 
modes  1  (panel  a)  and  2  (panel  6).  The  currents  are  computed  from  both  the 
beam  electrons  and  the  ambient  electrons  in  the  original  beam  region.  The 
dramatic  growth  in  all  four  cases  during  the  period  from  IScDpe'1  to  SOGOpe*1 
illustrates  the  common  origin  of  the  current  structure  in  the  processes 
associated  with  the  beam-plasma  interaction.  Although  the  current 
structure  has  significant  components  in  modes  1  through  7,  only  the  first 
two  can  satisfy  the  Cerenkov  resonance  condition  k||  =  co/vh  with  to  <  Qe  and 
vii  $  vb  and  thus  produce  whistler  waves.  Since  IJ||I2  »  IJj_l2  in  the 
simulations,  one  surmises  that  the  whistler  waves  are  generated 
predominantly  by  the  J||  current.  This  conclusion  is  further  strengthened  by 
the  single-particle  analysis  of  Mansfield  [1967]  which  shows  that  vjl  makes 
only  a  small  contribution  to  the  Cerenkov  radiation  level  if  vj.  <  v|| . 

The  initial-value  simulations  have  thus  demonstrated  that  whistler 
waves  are  produced  from  electron  beams  on  the  same  time  scale  as  the 
classical  beam-plasma  interaction  whose  dominant  feature  is  the 
production  of  intense  Langmuir  oscillations.  The  growth  of  these 
electrostatic  oscillations  is  accompanied  by  the  formation  of  a  current 
structure  which  acts  like  an  antenna  and  emits  whistler  waves  in  a 
coherent  manner.  Since  the  generation  efficiency  in  the  simulations  is  so 
much  higher  than  that  actually  observed,  however,  it  is  not  clear  that  this 
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Figure  11.  Time  histories  of  ljy(ky)l2  and  ljz(ky)l2  in  the  beam  region  for 
(a)  mode  1  and  (b)  mode  2. 


mechanism  is  actually  the  operative  one  in  the  beam  experiments  in  space. 
We  address  this  point  in  the  next  subsection. 

6.2.  Radiation  from  an  Injected  Beam 

We  now  use  the  electromagnetic  version  of  the  isolated-system 
simulation  code  to  reexamine  the  process  of  coherent  whistler  radiation 
from  electron  beams.  This  is  a  much  more  realistic  model  than  the  initial- 
value  calculation  since  the  beam  is  now  injected  at  a  single  point.  The 
configuration  is  similar  to  that  shown  in  Figure  1;  the  beam  energy  is  now 
higher  (Eb  =  16  keV),  and  the  grid  spacing  A  corresponds  to  about  40  cm. 
The  magnetic  field  strength  will  be  varied  but  is  always  such  that  the 
electron  cyclotron  frequency  Oe  is  less  than  the  ambient  plasma  frequency 
(Ope,  as  was  the  case  during  the  Spacelab  2  observations. 

The  most  striking  difference  between  the  present  continuous-injection 
simulations  and  the  previous  initial-value  results  is  the  sharply  reduced 
level  of  electrostatic  Langmuir  oscillations.  Instead  of  a  rapid  initial  growth 
followed  by  collapse  of  these  oscillations,  there  is  now  a  slower  build  up  of 
the  parallel  wave  energy  toward  a  quasi-static  level  on  a  time  scale  of  30-40 
oipb'1.  The  total  electrostatic  field  energy  is  now  only  a  few  tenths  of  a 
percent  of  the  energy  of  the  beam  electrons  injected  during  this  interval. 
The  initial-value  simulations  thus  appear  to  have  overestimated 
significantly  the  strength  of  the  beam-plasma  interaction. 

Figure  12  shows  contours  of  the  perturbed  Bz  field  at  times  copbt  =  40,  80, 
and  120  for  the  case  of  &e  /©pe  =  0.5.  The  current  sheet  formed  by  the  injected 
electron  beam  produces  an  antisymmetric  field  relative  to  the  beam  axis. 
Superimposed  on  this  dominant  field  are  periodic  modulations  with  a 
wavelength  of  ~  90A.  The  direction  of  propagation  relative  to  the  ambient 
magnetic  field  appears  to  decrease  as  a  function  of  time,  from  *  10°  at  the 
earliest  time  to  «  4*  at  the  latest  time.  The  perturbed  By  field  shows  a 
similar  wave  pattern.  Figure  13  shows  power  spectra  measured  at  a  point 
Ax  =  90A,  Ay  =  40A  relative  to  the  beam  injection  point  for  the  y  and 
z  components  of  the  magnetic  field  and  the  y  component  of  the  transverse 
(electromagnetic)  electric  field.  The  By  spectrum  shows  a  clear  peak  at 
co  ~  0.7  the  Ey  spectrum  has  a  similar  peak  but  with  an  intensity  about 
two  orders  of  magnitude  weaker.  We  thus  identify  the  wave  disturbances  as 


Figure  12.  Contours  of  the  Bz  magnetic  field  resulting  from  the 
injection  of  an  electron  beam  located  at  x/A  =  245  and 
y/A  =  64  into  an  ambient  plasma  with  Qe/c&pe  =  0.5.  Solid 
contours  indicate  positive  values  of  Bz,  while  dotted 
contours  indicate  negative  values. 
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Figure  13.  Power  spectra  for  the  By,  Bz,  and  ETy  fields  averaged  over 
the  interval  0  <  ©pbt  <  200  at  a  point  8x/A  =  90, 5y/A  =  40  away 
from  the  source  of  the  electron  beam.  The  ambient  density 
is  given  by  iV©pe  =  0*5- 


whistler  waves  propagating  with  a  phase  velocity  of  order  c/10.  The  Bz 
spectrum  shows  a  similar  spectral  intensity  in  the  whistler  region,  but  it 
produces  only  a  shoulder  on  the  much  larger  fields  associated  with  the 
current  sheet.  The  intense  peaks  at  frequencies  5  copi  are  associated  with 
the  ion  response  to  establish  current  closure  across  field  lines.  The  total 
whistler  wave  energy  is  now  about  5  x  10*5  of  the  electron  beam  energy 
injected  up  to  time  <opbt  =  80.  This  is  now  consistent  with  the  Spacelab  2 
observations.  Similar  results  are  obtained  for  a  second  simulation  with 
o)pe/Qe  -  0.3  except  that  the  angles  of  propagation  are  somewhat  larger, 
ranging  from  «  19*  down  to  «  8*  as  a  function  of  time. 

With  the  present  simulation  parameters  the  spacecraft  potential 
increases  slowly  with  time  (after  a  rapid  initial  jump  and  decay)  and 
reaches  a  level  of  *  40%  of  the  beam  energy  by  the  end  of  the  run  (b)pbt  =  200). 
As  a  result  of  this  charging,  the  ambient  electron  density  is  reduced  on  field 
lines  connected  to  the  spacecraft.  This  leads  to  the  formation  of  a  density 
depression  extending  out  from  the  spacecraft  whose  spatial  extent  along  the 
magnetic  field  increases  with  time  at  a  speed  of  «  vb/3.  The  average  density 
decrease  is  about  a  factor  of  two,  but  there  are  local  depressions  of  a  factor  of 
four  or  more.  The  establishment  of  this  density  trough  has  a  significant 
effect  on  the  propagation  of  the  whistler  waves.  For  0.5  <  oo/Qe  <  1,  whistler 
waves  can  be  trapped  in  a  density  trough  which  then  acts  as  a  duct  [see 
e.g.,  Smith  et  al.,  I960].  This  explains  the  observed  decrease  in  propagation 
angle  of  the  waves  as  a  function  of  time  as  the  trough  expands.  The  density 
depressions  produced  in  the  simulations  are  considerably  larger  than  the 
values  of  2-10%  necessary  for  trapping. 

7.  SPACECRAFT  POTENTIAL  CONTROL  VIA  INDIUM 
ION  BEAM  INJECTION 

A  recurring  difficulty  for  space  satellites  making  plasma 
observations  concerns  the  accurate  measurement  of  very  cold  plasma 
and/or  low  frequency  electric  fields  during  periods  of  high  electrostatic 
surface  potentials.  Although  scientific  spacecraft  have  been  equipped  with 
conductive  surfaces  to  eliminate  differential  charging  between  sunlit  and 
shadowed  surfaces,  the  fundamental  problem  of  overall  electrostatic 
charging  remains.  Olsen  [1982]  used  the  reduction  of  the  spacecraft 
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potential  during  a  transit  from  sunlight  into  eclipse  for  the  identification  of 
a  cold  plasma  population  previously  undetected  due  to  the  masking  effect  of 
the  positive  surface  potential. 

The  electrostatic  potential  of  the  spacecraft  is  a  consequence  of  the 
current  equilibrium  between  charged  particles  impacting  on  the  surface 
and  charges  leaving  the  surface,  such  as  photo-electrons,  sputter-induced 
charged  secondary  particles,  and  artificially  ejected  charges  [Pedersen  et 
al.,  1987].  The  extreme  values  of  the  electrostatic  potential  may  range  from 
less  than  -1000  Volts  up  to  +1000  Volts  [ Garrett ,  19811.  To  maintain  the 
spacecraft  potential  near  the  ambient  plasma  potential,  an  ion  emission 
system  that  is  small  in  size  and  resources  is  needed.  A  major  constraint  on 
such  a  system  comes  from  the  fact  that  its  spacecraft  resource 
requirements  remain  within  a  reasonable  ratio  with  other  instruments  on 
the  spacecraft  and  it  should  also  be  "quiet"  in  the  sense  of  not  interfering 
with  other  instruments  onboard.  Keeping  these  requirements  in  mind  it 
was  decided  to  develop  new  systems  based  on  the  use  of  liquid  indium 
[Schmidt  et  al .,  1989]. 

Included  in  the  instrumentation  package  for  the  upcoming  Cluster 
and  Geotail  satellites  will  be  a  small  compact  ion  emitter  that  contains  a 
tiny  amount  of  indium  carried  in  a  heater  controlled  containment  device. 
Once  the  indium  is  heated  to  well  above  its  melting  point,  the  liquid  metal 
flows  due  to  surface  tension  along  the  shaft  of  a  small  needle.  At  its  top  the 
indium  is  ionized  by  field  emission  and  emitted.  The  typical  beam  drift 
kinetic  energy  is  5  to  6  keV. 

The  main  concern  of  the  research  we  have  carried  out  has  been  to 
examine  whether  the  emitted  indium  ion  beam  can  be  unstable  to 
electrostatic  plasma  instabilities  for  the  various  emitted  beam  and 
background  conditions.  For  this  study  we  restrict  ourselves  to  a  standard 
size  spacecraft  with  conductive  surfaces  and  to  orbits  outside  the 
plasmasphere.  Eclipse  conditions  are  ignored  and  therefore  the  sunlit 
spacecraft  will  charge  positive  due  to  the  dominance  of  sun-induced  photo¬ 
electrons.  The  goal  of  the  linear  theory  study  is  to  identify  any  unstable 
wave  modes  that  can  occur  such  that  effects  that  may  result  from  these 
instabilities  can  be  distinguished  from  other  plasma  physical  processes. 
The  detailed  parameter  search  quantifies  under  what  conditions  the 
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indium  ion  beam  can  be  emitted  in  sucli  a  way  that  an  instability  can  be 
turned  on  or  off  during  the  mission. 


8.  OPERATING  PRINCIPLES  OF  THE  LIQUID  METAL  ION 
SOURCE  (LMIS) 

The  ion  emitter  [Riidenauer  et  al.,  19881  is  a  "solid-needle"  type  liquid 
metal  ion  source  using  indium  as  charged  material  [Dixon  and  von  Engle, 
1980].  Indium  has  an  atomic  number  of  49  and  its  mass  is  114.82.  A 
tungsten  needle  with  a  tip  radius  of  typically  2  to  15  p  is  mounted  in  the 
indium  reservoir  which  is  attached  to  a  heater.  If  the  needle  is  well  wetted 
by  the  metal,  then  the  electrostatic  stress  at  the  needle's  tip  pulls  the  liquid 
indium  towards  the  extractor  electrode.  This  stress  is  counteracted  by  the 
surface  tension.  The  resulting  equilibrium  is  characterized  by  a  stable  field 
emission  from  a  liquefied  cone  (Taylor  cone)  with  a  tip  radius  of  1  to  5  nm. 
The  evaporation  and  resulting  ionization  leads  to  the  emission  of  positively 
charged  ions  from  the  Taylor  cone.  Typical  current  densities  at  the  cone 
apex  are  of  the  order  of  106cm‘2sr*1A  for  a  beam  energy  of  10  keV.  Ions 
leaving  the  surface  are  constantly  replenished  by  hydrodynamic  flow  of 
liquid  indium  from  the  reservoir  to  the  tip  of  the  needle. 

The  typical  operating  characteristics  of  the  LMIS  require  an  extraction 
high  voltage  of  5  kV  applied  between  extraction  electrode  and 
needle/reservoir.  An  emission  current  of  10  pA  requires  a  supply  current  of 
about  20  pA,  mostly  due  to  focusing  losses.  The  heater  power  required  to 
maintain  the  indium  reservoir  well  above  the  melting  temperature  is  about 
0.4  W.  The  beam  divergence  at  the  orifice  is  better  than  15  (half-width),  and 
the  beam  temperature  AE/E  is  better  than  10‘3  at  1  pA,  with  a  weak 
dependence  on  the  applied  extraction  voltage.  The  beam  temperature 
depends  on  the  emission  current,  which  can  be  set  by  adjusting  a  current 
limiter  inside  the  high  voltage  power  supply.  For  small  emission  currents, 
around  1  pA,  the  measured  beam  temperature  is  of  the  order  of  40  eV. 
Increasing  the  emission  current  up  to  10  pA  results  in  an  increasing  beam 
temperature  (600  eV)  that  is  mostly  due  to  the  appearance  of  a  high  energy 
tail. 

For  the  linear  theory  study,  we  take  the  In+  beam  drift  energy  to  be  5 
keV,  which  corresponds  to  a  drift  speed  of  about  91.3  km/s.  The 
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temperature  of  the  beam,  as  explained,  is  variable  depending  on  the 
emission  current  setting  and  can  range  from  40  to  600  eV;  the  beam 
thermal  speed  will  be  a  parameter  that  will  be  varied  during  the  linear 
theory  study. 

9.  LINEAR  THEORY 

A  linear  theory  parameter  search  has  been  carried  out  to  investigate 
whether  particular  combinations  of  ambient  plasma  parameters  and 
emission  modes  of  the  indium  ion  beam  exist  such  that  electrostatic  plasma 
instabilities  may  be  excited.  We  consider  a  three  species  plasma  with  (1)  an 
indium  ion  beam,  (2)  background  electrons,  and  (3)  background  ions.  The 
background  plasma  parameters  vary  depending  on  where  the  satellite  is  at 
the  time  of  beam  emission.  For  our  study  we  considered  four  different 
background  regions:  the  solar  wind,  the  lobe,  the  plasma  sheet  boundary 
layer,  and  the  central  plasma  sheet.  Typical  plasma  parameters  for  each  of 
these  regions  used  in  this  study  are  found  in  Table  2,  and  are  based  on 
previous  observations  made  in  each  region  [Frank,  1985;  Lui,  1987; 
Rosenberg  and  Whipple,  1990]. 


TABLE  2.  Typical  observed  parameters  for  each  of  the  four 
regions  under  consideration 


Region 

Ti/Te 

COpe^le 

Te(eV) 

ne(cm‘3> 

CPS 

5 

5-20 

>1000 

1 

PSBL 

2 

5 

-500 

0.1 

Lobe 

2 

1 

<100 

0.01 

Solar  Wind 

<2 

55-200 

10 

1-10 

To  understand  the  basic  physics  of  unstable  wave  modes  possible  for 
this  type  of  configuration,  we  first  consider  a  cold  plasma  with  magnetized 
electrons  and  unmagnetized  ions.  Following  Mikhailovskii  [1975],  the 
dispersion  relation  can  be  written  as 
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(Qpe2  sin29  (Ope2  cos2  9  _  (Opj2  (0pb2 
CO2  -  fle2  (O2  (02  (ft)  -  kU)2 

where  (Ope,  (0pi,  and  ©pb  are  the  plasma  frequencies  of  the  background 
electrons,  background  ions  and  beam  ions,  respectively,  Qe  is  the  electron 
gyrofrequency,  0  is  the  angle  between  the  wavenumber  k  and  the  external 
magnetic  field,  and  l)  is  the  indium  beam  drift  velocity.  We  assume  that 
kand  C?  are  aligned.  The  frequency  (0  is  complex  with  co  =  (0r  +  iy. 

An  unstable  interaction  between  the  beam  mode  and  the  lower  hybrid 
mode  occurs  at  very  oblique  wave  propagation  (0  -  90*)  when  ©r  =  coih  ~  kU, 
where  the  lower  hybrid  frequency  is  given  by 


(Oih  =  ■ . . -  (3) 

V 1  +  (Ope2/0e2 

in  the  limit  of  unmagnetized  ions,  i.e.,  (opi  »  Qi-  The  growth  rate  for  this 
instability  is  given  by 


_ mhay3 _ 

[  1  +  (mi/me)cos20] 1/3 


(4) 


where  a  =  nbrni/nimb  «  1,  nb  and  mb  are  the  beam  density  and  mass,  ni 
and  mi  are  the  background  ion  density  and  mass,  and  me  is  the 
background  electron  mass.  The  wave  growth  is  due  to  a  fluid-type 
instability  driven  by  the  net  drift  between  the  beam  and  background  plasma, 
and  it  is  interesting  to  note  from  equation  (4)  that  the  growth  rate 
maximizes  when  0  =  90*. 

If  we  consider  the  case  where  the  indium  beam  is  not  strictly  cold, 
but  is  approximated  by  a  Maxwellian  velocity  distribution,  a  beam  driven 
lower  hybrid  instability  with  real  frequency  given  by  equation  (3)  is  still 
present,  but  the  growth  rate  is  modified  and  can  be  expressed  as 


/  U  \  2  arnih 
\vtb/  1  +  (mi/me)cos20 


(5) 


where  vtb  is  the  beam  thermal  speed.  As  in  the  cold  plasma  case,  the 
instability  maximizes  for  transverse  wave  propagation  (0  =  90s),  but  for  a 
warm  beam  the  growth  rate  is  sensitive  to  the  beam  drift  to  thermal  speed 
ratio  U/vtb,  such  that  as  the  beam  temperature  increases  for  a  given  drift 
speed,  the  growth  rate  will  decrease. 
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To  study  the  parametric  dependence  of  this  instability  for  the 
different  background  regions  and  beam  emissLxi  modes,  we  have  solved  the 
fully  magnetized  warm  electrostatic  dispersion  relation  numerically.  The 
form  for  the  general  electrostatic  dispersion  relation  can  be  found  in  Harris 
[1959]  or  Stix  [1962].  An  example  of  the  dispersion  relation  for  the  indium 
ion  beam  ejected  while  the  satellite  is  in  the  plasma  sheet  boundary  layer  is 
shown  in  Figure  14.  The  beam  drift  speed  to  thermal  speed  ratio  used  for 
this  figure  is  U/vtb  =  17  and  the  beam  density  is  nb/ne  =  0.5.  Plotted  is 
frequency  versus  wavenumber  with  the  real  frequency  shown  by  the  solid 
line  and  the  growth  rate  shown  by  the  dashed  line.  In  Figure  14,  the  wave 
propagation  angle  is  varied  and  the  results  confirm  that  found  from 
equations  (4)  and  (5),  whereby  the  growth  rate  maximizes  for  purely 
transverse  wave  propagation  (0  =  90*).  Keep  in  mind  that  the  beam  is 
directed  parallel  to  k ,  which  means  that  maximum  wave  growth  occurs  for 
the  beam  being  emitted  90°  to  the  external  magnetic  field.  The  real 
frequency  at  maximum  wave  growth  is  less  than  the  lower  hybrid 
frequency  given  in  equation  (3),  due  to  warm  plasma  effects  of  the 
background  ions. 

Now  maintaining  transverse  wave  propagation  (0  =  90°)  to  find 
maximum  wave  growth,  we  plot  in  Figures  15,  16,  17,  and  18  the  dispersion 
relation  for  beam  emission  in  the  central  plasma  sheet,  plasma  sheet 
boundary  layer,  solar  wind,  and  lobe,  respectively  (the  background 
parameters  are  listed  in  Table  2).  The  format  for  these  figures  is  the  same 
as  in  Figure  14,  and  the  beam  density  is  the  same;  however,  the  ratio  U/vtb 
is  varied  (keeping  U  fixed  corresponding  to  a  beam  drift  energy  of  5  keV)  for 
each  region  to  find  the  instability  threshold.  As  expected  from  equation  (5) 
the  growth  rate  decreases  as  U/vtb  decreases,  which  corresponds  to  an 
increase  in  the  beam  temperature.  From  Figures  15  through  18  it  can  be 
seen  that  the  central  plasma  sheet  has  the  lowest  threshold  for  the  beam 
instability  as  a  function  of  beam  temperature,  which  means  that  the 
indium  beam  can  be  fairly  warm  and  still  drive  an  instability  there.  The 
lobe,  on  the  other  hand,  has  the  highest  threshold  whereby  the  beam  must 
be  relatively  cool  (low  emission  current)  for  an  instability  to  exist. 
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Plasma  Sheet  Boundary  Layer 


Figure  14.  This  figure  shows  the  dispersion  of  low  frequency  waves 
found  from  numerical  solutions  of  the  general  warm 
plasma  dispersion  relation  for  an  indium  ion  beam  emitted 
while  the  satellite  is  in  the  plasma  sheet  boundary  layer. 
Plotted  is  frequency  (normalized  to  the  background  ion 
gyrofrequency)  versus  wavenumber  (times  the  electron 
Debye  length).  The  solid  line  corresponds  to  the  real 
frequency  and  the  growth  rates  are  shown  by  the  dashed 
line  for  three  different  wave  propagation  angles.  Note  the 
growth  rate  maximizes  for  transverse  wave  propagation. 
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Central  Plasma  Sheet 
Dispersion  Relation 


kAe 


Figure  15.  The  wave  dispersion  relation  when  the  indium  beam  is 
emitted  in  the  central  plasma  sheet.  The  format  is  the  same 
as  Figure  14,  but  here  the  ratio  U/vtb  is  varied  and  the 
lowest  value  where  there  is  still  an  appreciable  instability  is 
shown.  The  background  plasma  parameters  are  found  in 
Table  2.  In  this  case,  if  U/vtb  falls  below  6,  no  instability  is 
possible. 
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Lobe  Dispersion  Relation 


14.1,  14.5,  15.3 


Figure  18.  Dispersion  relation  for  beam  emission  in  the  lobe.  The 
format  is  the  same  as  in  Figure  15,  with  the  background 
parameters  found  from  Table  2.  Note  that  this  region  has 
the  highest  threshold  for  instability  compared  to  the  others. 


10.  SUMMARY 


An  indium  ion  beam  device  will  be  placed  on  future  satellites  such  as 
Geotail  and  Cluster  to  control  the  spacecraft  potential.  A  study  has  been 
carried  out  to  see  if  the  emitted  indium  ion  beam  can  be  unstable  when  the 
satellite  is  in  various  regions  of  the  Earth's  magnetosphere  and  solar  wind. 
Linear  theory  for  electrostatic  waves  shows  that  an  ion  beam  driven  lower 
hybrid  instability  can  be  excited  with  maximum  growth  at  wave 
propagation  (and  beam  emission)  transverse  to  the  external  magnetic  field. 

Numerical  results  from  the  linear  theory  parameter  search  show 
that  for  the  expected  indium  beam  conditions  and  typical  parameters  for 
the  background  plasma  that  the  central  plasma  sheet  will  have  the  lowest 
threshold  for  instability,  while  the  lobe  will  have  the  highest.  Since  the 
emission  current,  which  controls  the  beam  thermal  speed,  can  be  varied 
during  the  satellite  mission,  the  linear  theory  results  can  be  tested  directly 
by  changing  the  current  and  using  observed  wave  data  to  see  whether  an 
instability  in  the  right  frequency  range  is  excited  and  turns  on  or  off 
according  to  the  stability  criteria.  Also,  depending  on  the  orientation  of  the 
beam  emitter  with  the  external  magnetic  field,  the  linear  theory  conclusion 
that  maximum  wave  growth  will  occur  for  transverse  emission  can  be 
tested. 
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